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Abstract 

We present a new, completely revised calculation of the 
muon anomalous magnetic moment, = [g^ — 2)/2, comparing 
it with the more recent experimental determination of this quan- 
tity; this furnishes an important test of theories of strong, weak and 
electromagnetic interactions. These theoretical and experimental 
determinations give the very precise numbers, 

r 116 591 806 ± 50 ± 10 (rad.) ± 30 (£ x I) [Th., no r] 
lO" xa^, = \ 116 591 889 ± 49 ± 10 (rad.) ± 30 x I) [Theory, t] 
I 116 592 080 ± 60 [Experiment]. 

In the theoretical evaluations, the first quantity does not, and the 
second one does, use information from t decay. The first errors for 
the theoretical evaluations include statistical plus systematic errors; 
the other ones are the estimated errors due to incomplete treatment 
of radiative corrections and the estimated error in the light-by-light 
scattering contribution. We thus have a significant mismatch be- 
tween theory and experiment. We also use part of the theoreti- 
cal calculations to give a precise evaluation of the electromagnetic 
coupling on the Z, qq.e.d .{Ml), of the masses and widths of the 
(charged and neutral) rho resonances, of the scattering length and 
effective range for the P wave in tttt scattering, and of the quadratic 
radius and second coefficient of the pion form factor. 
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1. Introduction 



The anomalous magnetic moments of electrons and muons provide one of the more im- 
pressive tests of the standard model of strong, weak and electromagnetic interactions. The 
electron anomaly, a^, receives only marginal contributions from weak and strong interac- 
tions, being dominated by electromagnetic (QED) radiative corrections. The agreement 
between theory and experiment is such that, if we turn it around, the experimental value 
of Oe provides the more precise determination of the fine structure constant, t^l a. 

For the nmon magnetic moment, the recent precise measurements!'^! Qf the muon 
anomaly, a^u, constitute one of the more impressive tests not only of electroweak inter- 
actions, but of strong interactions as well. After these measurements, the world average 
value for = {g^jt — 2)/2 is 

10^^ X = 116 592 080 ± 60 [Experiment]. (1.1) 

Considering that theory also gives the value of the magnetic moment itself (and not only 
the anomaly) the agreement of theory and experiment that we will describe represents a 
precision of one or two parts in a billion. Nevertheless, and as we will see, there remains 
a discrepancy at the level of 2.3 a to 3.3 cr. 

The electromagnetic and weak contributions to have been calculated with great 
accuracy:'^! 

10^^ X a^(QED) = 116 584 719 ± 1.8, 

10^^ X a^(Weak) = 152 ± 3. ^^'"^^ 

Combining this with (1.1), we find the experimentai number for the hadronic contributions 
to a^, 

10^^ X a^(Hadr.) = 7209 ± 60 [Experiment]. (1.3) 

The evaluation of this quantity, a^(Hadr.), from theory will be the main subject of the 
present note; we will find 

1 0il X « (yr. . _ r 6 935 ± 50 ± 10 (rad.) ± 30 (£ x i) [No r] 

iu x«Mliiaar.;- <j^^p^g^^g^^Q^^^^^^^gQ^^^^^ [With r]. ^ > 

The numbers above depend on whether or not one includes information on r decay, which 
is probably the more reliable result as it is the one that incorporates more information: see 

our text below for details. In the errors, "rad" and "£ x ^" refer, respectively, to estimated 
errors due to uncalculated radiative corrections and light-by-light scattering contributions. 
The results in (1.4) takes into account all the more recent e+e^ annihilations data, and 
CTT scattering data. At the end of the present article we will comment on the degree 
of agreement of theory and experiment, compare our results with those of other recent 
calculations, and discuss the possible reasons for the discrepancy between (1.3) and (1.4). 

We will also give a summary of the results that some of the calculations imply 
for the electromagnetic coupling on the Z, aQ.E.D.(^z)' masses and widths of the 

(charged and neutral) rho resonances, for the scattering length and effective radius for the 
P wave in pion-pion scattering, and for the quadratic radius and second coefficient of the 
pion form factor. 
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Figure 1 . The order o? 
hadronic contributions to the 
muon magnetic moment. The 
blob represents an arbitrary 
hadronic state. The wavy 
hnes are photons. 



2. The hadronic contributions to a^^. I: the 0{o?) piece 



To order o? ^ the contributions to a^(Hadr.) can be represented by the diagram shown in 
Fig. 1. As has been known for a long time, they can be written in terms of the cross 
section for e+e" annihilation into hadrons, as follows. We first write 



/"OO 

a(f)(Hadr.) =1277 / As K{s)\m.n{s) 



K{s)=^K{s); K{s)= dx^—^ -i^. 

37r^s Jo + (1 — x)s/m-^ 

Here 77 is the hadronic part of the photon vacuum polarization function. Then, we can 
express Im77 in terms of the ratio of (lowest order) cross sections for e^e~ annihilation 
into hadrons over annihilation into muons: 

(t(°' (e+e^ ^ hadrons) rni / -i- - -i- -x 47ra^ 



a(f)(Hadr.) = / dsii:(s)i?(s). (2.1b) 

Therefore, the situation is, in principle, simple: we take the experimental cross 
hadrons, insert it into (2.1b) and the value for a)j (Hadr.) will follow. 
In practice, however, things are more complex. The experimental numbers for the cross 
section e'^e~ into hadrons present (relatively) large errors and we have, therefore, interest 
to supplement these with more precise theoretical formulas whenever possible. As a matter 
of fact, by so doing we are able to diminish the error in the theoretical calculation almost 
by a factor of two. 
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2.1 The low energy region, s < 0.8 GeV^ 

In the region below s = 0.8 GeV^ that wc may caU "rho region'' , experimental data have 
in fact improved substantially in the last years due to the more recent measurements, 
especially at Novosibirsk, t*'^' Here the important contributions are those of the omega 
resonance (for which experimental data and the Goimnaris-Sakurai method may be used; 
see below and ref. 6) and the two-pion contribution, which is the one that may be made 
more precise using theory, and also the one that we will discuss in more detail because 
it provides the bulk of a^(Hadr.). Here one profits from the fact that the two pion 
contribution can be expressed in terms of the pion form factor, F„ , 

1 / 4rv,2 \ 

Iran,^is) = — [l-^j \F^is)f, (2.2) 

and may be determined from fits to data on e+e^ — > 2tt and, using analyticity, also 
data from ttc — > Tre scattering, i.e., at spacelike s. What is more, we may use (with due 
caution; see below) data on t decay, — > vtt^tt^^ related to F„ by isospin invariance. 

We will for the moment work in the approximation of neglecting electroweak 
corrections to F^r; we will discuss this further in Subsect. 2.2. In this approximation, the 
properties of F^^ (s) that allow us an improved calculation are the following: 

(i) F^(s) is an analytic function of s, with a cut from 4m^ to infinity. 

(ii) On the cut, the phase of F.;r(s) is, because of unitarity, identical to that of the P- 
wave in tttt scattering, 5i(s), and this equality holds until the opening of the inelastic 
threshold at s = to (Fermi- Watson final state interaction theorem) . 

(iii) For large s, Ft^{s) ~ 1/s. Actually, one knows the coefficient of this behaviour, but 
we will not need it here. 

(iv) F^(0) = 1. 

The inelastic threshold occurs, rigorously speaking, at s = 16rn^. However, it 
is an experimental fact that inelasticity is negligible until the quasi-two body channels 
un, oiTT . . . are open. In practice one can take 

to ^ 1 GeV^ 

and fix the best value for to empirically. It will he to — 1-1 GeV'^, and the dependence of 
our results on to is very slight. 

The properties (i-iv) can be taken into account with the Omnes-Musklielishvili 
method.^ We construct a function J(s) with the proper phase and asymptotic behaviour 
by defining 



/ „ \ [l-Sl(to)/n]to/s / 



exp 




1 



More details about the solution of the Omnes-Muskhelishvili equations can be found in the 
text of N. I. Muskhelishvili, Singular Integral Equations, Nordhoof, 1958, and, applied to our 
present case, in ref. 6. 
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We have written the dispersion relation with one subtraction to ensure that J(0) = 1. 
The singular integral is understood to be calculated replacing s — > s + ie, e > 0, e — > 0. 
We then define the function G by 

F^(s) = G(s)J(s), (2.3b) 

and it follows from properties (i-ii) that G{s) is analytic with only the exception of a cut 
from to to infinity, as we have already extracted the correct phase below s = to- 

We can apply the effective range theory to parametrize the phase 5l. According 
to this, the function 

^(t) cot 5l{t), fc=V^^M (2.4a) 

is analytic in the variable t except for two cuts: a cut from — oo to 0, and a cut from 
t — to to +00. To profit from the analyticity properties of V' we will make a conformal 
transformation. We define 

w = —= 

Vt + Vt^ 

When t runs the cuts, w goes around the unit circle. We may therefore expand ^ in a 
power series convergent inside the unit disc. In fact, because we know that the P wave 
resonates (which implies a zero of xp) it is convenient to expand not i/j itself, but the ratio 
ip{t)/{mp — t) = ip{t). Here rUp is the mass of the rho resonance; so we write, 

^{t) = {ml - t)i>{t) = {ml - t) {bo + biw + ■■■}. (2.4b) 
The P-wave, 7 = 1 tttt scattering length, aj, is related to tp by 

^ mjr?/'(4mj) ' 

experimentally, a\ ~ (0.038 ± 0.003)7n^^, a condition that may be incorporated into the 
fit. Note, however, that we do not assume the values of mp, Fp. We only require that ip 
has a zero, and will let the fits fix its location and residue. It turns out that, to reproduce 
the width and scattering length, and to fit the pion form factor as well, only two bo, bi 
are needed in (2.4b). 

We now turn to the function G{s). This function is analytic except for a cut from 
s = to to +00. The conformal transformation 



h^/to — Vto' 

z = 



IVh + Vto - s 

maps this cut plane into the unit circle. So we may write the expansion, 

G{s) = l + ci{z + 1/3) + C2{z^ - 1/9) + • • • , (2.5) 

which takes into account the condition G(0) = 1 order by order. We will only need two 
terms in the expansion, so we have ci, C2 as free parameters. This means that, altogether, 
we have the five parameters, 

nip, bo, bi, ci, C2, 
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to fit 158 experimental points. 

One can then use the formulas just discussed and fit the experimental data on 
after we have taken into account the uj — p interference (which includes the a; — > 27r 
piece). This we do with the Gounnaris-Sakurai method. We write 

Ml 

1+CT- 



(s)=F^-(s)x M^=m^-\rj2 (2.6a) 

where F^^^^ is the form factor we would have in absence of a; — p interference and the 
(real) omega mass. We take the values of , from the Particle Data Tables and find 

|ct| = (18±1) X lO"'^, argcT = 12±3°. (2.6b) 

The fit improves when using the 2002 Novosibirsk data;W there is better agreement with 
the data of Barkov et al.,1'*^ in particular for the larger values of s. The parameters of the 
fit are, however, similar to what we found in ref. 6 using the old (1999) Novosibirsk data, 
as we shall see. It should be noted that, in these data, electromagnetic corrections have 
been extracted; we give the details of the procedure in Subsect. 2.2 below. 

The result of the contribution to the hadronic part of the muon anomaly is now, 
fitting the e+e" annihilation data,[^'^l 

10^^ X ajf^(s < 0.8 GeV^) =4 707 ±21; xVd-o.f. = 91/(114 - 7). (2.7) 

The fact that the x^/d.o.f. is substantially smaller than unity means that there is some 
room for displacement of the central value given in (2.7). 

The result in (2.7) number may be improved, which we will do in two steps. First 
of all, we remark that, because our formulas for F^(s) are valid for spacelike as well as 
timelike s, we can use information not only from e+e^ annihilation, I'*'''! hut also from 
7re scattering, I ''^ which gives F^(s) for negative s. We record the results of two recent 
evaluations, made using this method. We have our evaluation here, 



fe^e , Trel 



10" X a(f)(s < 0.8 GeV^) 



172 

4 750 ±19 [TY; fixed norm.]; ^ - 



d.o.f . 159 - 7 

134 



4 715 ±20 ±25 [TY; float, norm.]; , ^ 

d.o.r. 161 — 10 

(2.8a) 

The errors in the first number here do not include systematic errors in ctt scattering, which 
is the reason for the largish x^/d.o.f . In the second set, the first error is the statistical error, 
the second the systematic one. We take into account the systematic errors, both in e^e" 
data!'*! and in the en data'^l as in ref. 6, by including a factor 1 + e in the normalization 
of each set of data, and allowing e to float. We flnd 



^eTT [ref. 7] — (1-3 ± 0.2)%, Ce+e- [Akhmetshin, ref.4] — (—0.4 ± O.t^ . 

(2.8b) 

^e+e- [OLYA; Barkov, ref.4] — (0-5 ± 0.9)%; 
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50 
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t (GeV^) 

Figure 2. The fit to the pion form factor data in the timehke and 
spacehke regions; the more recent (2002) Novosibirsk data are in- 
cluded. 

the fit is depicted in Fig. 2. The only e that is not compatible with zero is CgTr- In fact, we 
left its value arbitrary; that the result we found is consistent with the value quoted in ref. 7 
itself, €e7r — 1%, is a nontrivial test of the quality of data in this reference (in particular, 
of their estimate of systematic errors) and of the consistency of our fitting procedure. 
Then, we have a recent result by Colangelo and collaborators, I®' 

[e+e-, 7re] : lO" x a^^\s < 0.8 GeV^) = 4 679 ± 30 [Colangelo, ref. 8]. (2.9) 

The evaluation of Colangelo and collaborators, Eq. (2.9), includes the four-pion cut, and 
imposes the P wave phase shift as given by chiral-dispersive evaluations. This last fea- 
ture, however, makes the Colangelo result vulnerable to possible defects of the chiral- 
dispersive evaluation, such as those discussed in ref. 10. 
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A way to improve the precision of (2.7) is to include information on r decay^^^l in 
the fit. We consider the weak vector current correlator, related to r decay, and expand it 
as 



n^,=ij d^xe'P--{0\TV;{x)V.mO) = {-p^g^. +p^p,) n"" (t) + p^p.n^ (t); 



Vi=2TTlmn^, t 
We can then write 



^1 = 12 



1 



(m^+ -m^o) 



21 



1 - 



(m^+ + m^o) 
i 



3/2 



|FM(i)|2. 



(2.10a) 



(2.10b) 



To compare with the experimentally measured quantity, which involves all of Imi7^, we 

have to neglect the scalar component 77'^, which is proportional to (m^ — m„)-^, and thus 
very small. Moreover, F^^^ only equals F„ in the limit of exact isospin invariance. We 
have also complications due to radiative corrections, that we will discuss in Subsect. 2.3. 

One can take into account isospin breaking effects, at least partially, by realizing 
that the mass and widths of the and can be different. One finds. 



[e ' e , Tre, rj 



10 



11 



a''i^\s <0.8GeV^) = < 



4 793 ± 17 (St.), 



d.o.f. 



4 798 ±17 (St.) ±25 (Sys.) 



283 
241-9 



; [TY] 



245 



[TY, Sys.] 



d.o.f. 244- 13' 

(2.11a) 

depending on whether or not one takes into account systematic normalization errors. We 
unify the normalization of the tau decay data taking into account the relevant branching 
ratios, as given in the Particle Data Tables.^ 

In the second set of numbers in (2.11), we have, as in (2.8), taken into account 
systematic errors by allowing floating normalization by a factor 1 + e of the various data 
sets. We find, 

eew [ref. 7] = (1-0 ± 0.2)% [1%], [Akhmetshin, ref.4] = (0.4 ± 0.5)% [0.6%], 

= (-1.4±0.5)% [0.7%], ee+e- [OLYA; Barkov, ref. 4] — 

(2.0 ±0.8)% [4%]. 

(2.11b) 

In square brackets we give the estimate of the normalization errors as given by the exper- 
imental groups themselves, except for r decay, where the number 0.7% is taken from the 
Particle Data Tables. The only case where the value we find for e exceeds the expectations 
is for T decay data, although the difference is very small, 0.7 ± 0.5%. We will discuss this 
again in Sect. 2.3. 

The best number in (2.11a), of course, is the second, 4 798 ± 31. The fact that 
the Y^/d.o.f. is still a bit larger than unity can be traced partly to the size of e,-, that we 



^ When quoting the PDT we refer to the 2002 edition, K. Hagiwara et al., Phys. Rev. D66, 
010001 (2002). 
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discuss in Subsect. 2.3, and partly to a discrepancy between the tau decay data of OPAL 
and of ALEPH and CLEO. This is seen very clearly if we give the individual values for 
the ratio of to number of experimental points of the various sets of data, which we do 
for the second fit in (2.11a), i.e., including systematic errors: 

eTT [NA7, ref. 7] : 42/45 

e+e- [ref. 4] : 108/113 
r decay [Aleph, ref. 11] : 19/21 (2.11c) 
r decay [Cleo, ref. 11] : 32/30 
T decay [Opal, ref. 11] : 40/31. 

A remarkable feature of (2.11c) is that it shows that including information from tau decay 
in the fit does not spoil the quality of the fit of the pure e+e" data. 

2.2. Radiative corrections for the e^e~ Ccise 

We next devote a few words to discuss electromagnetic radiative corrections, a subject of 
crucial importance given the precision of the more recent data, but not always very clear 
in the existing literature. 

There are, in fact, two separate questions here. First, we have the radiative 
corrections to the hadronic part of the photon vacuum polarization, which affect Fig. 1 
by adding photon corrections (depicted below, in Fig. 6). These give corrections to 
of order c? and will be considered later, in Sect. 3. Secondly, we have the matter of the 
radiative corrections that have to be taken into account when extracting the pion form 
factor from experimental data. These come about for the following reason: the form factor 
that verifies the analyticity and unitarity properties necessary to carry out our analysis in 
Subsect. 2.1 is defined by 

WM\V') =i2n)-^p-p')^F„i{p-p'f), (2.12) 

clcctrowcak int.^0 

i.e., only strong interactions are taken into account for the expectation value of the elec- 
tromagnetic current, J^. 

This quantity, F^, has thus to be extracted from the experimentally measured 
cross sections for e+e~ — > tt+tt", which include all sorts of radiative corrections. To first 
(relative) order in a, these are the following: I) Corrections to the e''"e~7 vertex, or photon 
radiation by the incoming e+ , . These are pure QED effects readily calculated and taken 
into account as a matter of course in experimental analyses. II) Vacuum polarization 
corrections to the photon propagator. These are known in terms of ImTT, and are also 
be subtracted easily. They are explicitly taken into account in the second paper in ref. 4. 
Ill) Corrections to the 7r+7r^7 vertex, or photon radiation by the outgoing pions (Fig. 3). 
We now say a few words about the last. 

Radiation of hard photons by the outgoing pions, as in Fig. 3B, is excluded by the 
experimental cuts applied when analyzing e+e^ tt+tt" scattering, which require the 
angle between the momenta of 7r+, vr^ to be close to 180°. One is thus left with the soft 
photon radiation and vertex correction shown in Fig. 3A. They can be calculated under 
the assumption that one can factorize the pion form factor and, given the actual values of 
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(A) (B) 

Figure 3. Diagrams subtracted for evaluating the pion form 
factor contribution in aj, {s < 0.8 GeV^), but included in the 
0{a^) contribution to a^(Hadr.). The blob represents the pion 
form factor, to zero order in electroweak interactions, as defined 
in Eq. (2.12). 



the experimental cuts applied to the momenta of the 7r+ , tt^ , this correction turns out to 
be very small. In this approximation, the corrections have been evaluated long ago and 
are, fortunately, explicitly given and extracted in the 2002 version of the Novosibirsk data 
(ref. 4). 

Prom this analysis it follows that, for the e'^e~ — > n'^ir~ case, one can subtract 
all corrections and really obtain F^, as defined in Eq. (2.12), from data with an error that 
is only of order o.-^. This is the quantity to which we can apply our theoretical analysis, 
as we did at the beginning of the present Section. 

A problem with the evaluation of the radiative corrections here, however, is that, 
as noted, one is using a model with elementary pions, in which the form factor is included 
by hand (factorized) . This may cause errors (for example, due to rescattering of the pions 
or dependence of the form factor on the off-shell pion mass), whose size (likely small, since 
the correction itself is small) we will estimate when discussing the case of tau decay in 
next Subsection. 

We should here note that the radiative corrections to the 7r"'"7r~7 vertex, or pho- 
ton radiation by the outgoing pions will have to be considered again when we consider 
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contributions to due to photon vacuum polarization of 0(a;"^). This we will do in 
Sect. 3. 

2.3. Comment on combining results using e'^e~ — > tttt and r — > VTVir, and 
on radiative corrections for r decay 

In some recent papers much ado is made about the difference in the pion form factor, 
in the region near and above the rho resonance, depending whether it its extracted from 
e+e^ TTTT or from r vh-k. In fact, this difference is expected. To begin with, |F(s)p 
grows almost by a factor 50 around the rho.^ Thus, even a small difference between the 
masses and widths of and will imply a large difference in the form factors. This 
matter was studied carefully in ref. 6, which discussion we summarize and update now. 

The values for the anomaly in the rho region that one obtains depending what 
one fits, is, with only statistical errors, 

10" X a(f ) (. < 0.8 GeV^) ^ [VZtf. \ f""" " '''' (2-13) 
^ ~ ^ (_ 4820 ±11 [ from r —i- z/TTTT only]. ^ ' 

For the tau, we have X^/d.o.f. = 87/(83 — 5). The two numbers are well outside each 
other error bars and it is, therefore, dangerous to combine them in a direct manner. What, 
however, one can do, is to fit simultaneously e+e" — > tttt and r — » z/tttt data allowing for 
different values of the parameters nip, b^, bi, (and therefore, also different widths'*) but 
with the same Omnes-Muskhelishvili function G{s) in Eq. (2.3b). As discussed in ref. 6, 
we expect isospin breaking effects to be small for G(.s) since its imaginary part is different 
from zero only for s > 1.1 GeV^. 

Of course, when calculating ap(s < 0.8 GeV^) one uses the parameters rUp, bo, bi 
determined from the fit to e+e^ tttt, even if we use tau decay data to help fix G: we 
should perhaps emphasize that the result reported in (2.11) is not an average of e+e" and 
T results, but an evaluation of e+e^ tttt, using information on G{s) from tau decay. 

This use of different masses and widths for p"^ , p^ is not sufficient to remove the 
discrepancies between the form factors obtained from e'^e~ and r decay. A reason for at 
least part of the remaining discrepancy is that the experimental number given in ref. 11 
for r decay includes also the radiative decay; that is to say, one does not measure the 
quantity r{T — > i/tttt) but, in fact, r{T — > vmr) + r{T — > i/tttt + 7) + higher orders. So we 
should discuss radiative corrections also for tau decay. 

First of all, we must make more explicit what we used in our fits. We take the 
experimental numbers for the decay r — > z^tttt including an eventual photon. The corre- 
sponding width we denote by -Texp. (t i^ir^ TT^\'~f)) , and we divide this by the experimental 
decay rate for r — > Vr^Ve^ plus eventual gammas, I^exp.(T i^rS-Veil))- However, and as 
in the case of e^e~ annihilation, we would like to have the quantity F^^^ with electroweak 

^ There is also a significant difference between data from e+e^ and tau decay for 0.8 GeV^ < 
s < 1.1 GeV^ which, however, affects very little the result for ap.{s < 0.8 GeV^). 

* We take, however, equal P wave scattering lengths for ■K'^Tr~ and 7r°7r^. We have checked that 
the influence of this is negligible. 
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V 




n n 



Figure 4. Radiative corrections to r — > vktt decay. 



interactions set to zero. Therefore, we should correct the ratio 

rexp.(r^^^7r^7r°(7)) 

to obtain the quantity 

r(0)(T ^ lyreve) 

where the F^^^ are evaluated to lowest order in electroweak interactions. 

The corrections necessary to do this are as follows. We first have radiative cor- 
rections to the leptonic width, that give 

rexp.(r ^ l^rei^eil)) = (1 + ^e) T^") (r ^ i/7r±^"), (^e - (f " TT^) ^ ~ -0.004. 

This correction was incorporated in the analysis of Subsect. 2.1. 

Then we have the corrections to the liadronic width. We here have corrections 
similar to those in Fig. 3B; and diagrams similar to those in Fig. 3A: see Fig. 4. We start 
with the first, that is to say, the sum of loop corrections (plus radiation of a soft photon). 
As far as we know, this has not been calculated exactly. However, we may expect this 
to be dominated by the short distance piece, since this last contains a large logarithm, 
log Mz/rriT- This gives the correction 

1 + log ~ 1.019 (2.14) 
TT rrir 

as has been known from a long time.I^-^1 This piece we have also extracted in our analysis 
above. After so doing, there still remains a piece of the loop correction and a correction 
due to soft photon radiation. Since we have extracted the large log Mz/rrir piece, we 
expect this to be comparable to the like piece in e'^e~ — > tttt case, and thus small. 
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It remains to correct for the ratio 

r{T — > VTrn) + r{T — > VETT + 7) 

r(o)(T^i/7r7r) 



^+5.^ ' (2-15) 



with 7 a hard photon, i.e., a photon with energy larger than a given Eq. This, again, is 
not known, but one can approximate it by the infrared logarithmic piece, which gives an 
average correction 

~ — log ~ 0.008, (2.16) 

the last for 20 MeV < -Bq < 80 MeV, and we have neglected terms of 0{mTr/mT)- M^^^ is 
the average invariant mass of the two pions, that we take equal to rUp. Since this correction 
plus the remainder (after extracting the logarithm) of the loop correction are not known 
exactly, we have not included the correction (2.16) in our evaluations above: (2.16) will 
be part of the normalization factor e,-. 

When we allowed variations of the normalization for tan decay data,!-'^^! multi- 
plying their numbers by 1 + e, and letting e float, we found the excellent fit reported in 
(2.11) with = —1.4%. It unfortunately is not possible to understand all of this as due 
only to the neglect of radiative corrections; we have verified that, including 5^ as given 
in (2.16) in the fit only changes Et- to —1.2%. We have to admit that there is a residual 
discrepancy with the value for the error in normalization, |e^| — 0.7%, given in the Particle 
Data Tables; but, the difference is small and, indeed, we get a fit to all experimental data 
with a ^ per experimental point which is essentially unity. 

In spite of this, it is clear that here we have not fully determined radiative cor- 
rections; only the logarithms log Mz, logE'o are exact in Eqs. (2.14, 16). Moreover, the 
function G{s) will not be exactly invariant under isospin and, finally, a correction due to 
the fact that m„ / nid, although likely very small, also exists beyond the rho. We must 
thus conclude that our partial ignorance of isospin violations, very likely dominated by 
radiative corrections, implies a possible shift of the central value of the anomaly. A conser- 
vative estimate for these effects would be the difference between the two values obtained 
leaving the r decay normalization fixed, and the same allowing it to float. This gives the 
number 

Aad [lO^^ X a[f)(s < 0.8 GeV^)] ~ 10. (2.17) 

We will accept the same error for the process e+e" — > tttt, although it is probably smaller 
here. 



2.4. The region with s > 0.8 GeV^ 



For low energies (say, s < 2 GeV'^), and near quark thresholds, there is no alternative 
to using experimental data.'^^' Between 0.8 GeV^ and 1.2 GeV^ we fit the experimental 
dataW for the tt+tt^ channel; its contribution is (229 ± 3 ± 3) x 10~^^. For other final 
states [KK, 37r, 47r, . . . ) we use the e+e^ data of refs. 5 and 13, with the methods of 
ref. 6. For higher energies, s > 2 GeV^, and away from quark thresholds, we can use QCD 
formulas, [^^1 taking into account the more recent values of the masses of the quarks as well 
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Figure 5. A typical dia- 
gram for the hadronic light 

by light contributions to the 
muon magnetic moment. 



as the strong coupling, as- For the QCD calculations we take the following approximation: 
for n / massless quark flavours, with charges Qf, we write 



■ \ " ^2 



a 



Qj<l + — + (1.986 - 0.115n/) — 



+ 



6.64 - 1.20n/ - 0.005n^ - 1.240-=^—^ 



To this one adds mass and nonperturbative corrections. We take into account the 0{m^) 
effect for s, c, b quarks with running masses TOi(s), which correct R'^^^ by the amount 

-3 Yl Qlfh^iis)U + 28^ + {294.8 -l2.3nf)(^y\s-\ 



i=s^c,b 

The details may be seen in refs. 6, 15. 

Adding all the contributions, one has 

.(2 



W X al^'i-s > 0.8 GeV^) = 2 134 ± 35. 



(2.19) 



This number is almost the same that may be found in ref. 6, using the 1999 Novosibirsk 
data; including the new data has very little influence in this region. 

3. The hadronic contributions to a^. II: the O(q^) pieces 



A contribution in a class by itself is the hadronic light by light one, that we label "£ x 
So we split 



a(Hadronic, 0{a^)) = a(-One blob' hadronic, 0{a^)) + a{£ x £). 



(3.1) 
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We will start by considering the last, given diagrammatically by graphs like that of 
Fig. 5. This can be evaluated only using models. One can make a chiral model calculation 
(essentially, replacing the hadronic blob by the lightest hadronic state, a 7r°), or one can 
use a constituent quark model in which we replace the blob in Fig. 5 by a quark loop. 
For the chiral model calculation we have to introduce a cut-off, since the 7r° contribution 
diverges for large virtuality of the photon lines. The result depends on the cut-off (for the 
chiral calculation) or on the constituent mass chosen for the quarks. After the correction 
of a sign error in the evaluations of ref. 16 (see ref. 17) we find 

10" X a{e X ^) = 86 ± 25 [Chiral calculation]. (3.2a) 

Earlier calculations with the model, using VMD to cure its divergence, gave (HKS, 
ref. 16) 

10" X a{e X ^) = 52 ± 20 [tt^ pole (HKS)]. (3.2b) 

One could also take the estimate of the tt*^ pole from Hayakawa, Kinoshita and 
Sandal^^l and add the constituent quark loop, in which case we get 

lO" X a(£ X ^) = 98 ± 22 [Quark const. model+ pion pole]. (3.2c) 

One expects the Tr'^-dominated calculation to be valid for small values of the virtual photon 
momenta, and the constituent model to hold for large values of the same. Thus, almost half 
of the contribution to a{i x i) in the chiral calciilation comes from a region of momenta 
above 0.5 GeV, where the chiral perturbation theory starts to fail, while for this range of 
energies, and at least for the imaginary part of (diagonal) light by light scattering, the 
quark model reproduces reasonably well the experimental data, as measured in photon- 
photon scattering. In view of this, we will take here the figure 

lO" X a(^ X ^) = 92 ± 30, (3.3) 

but will refrain from combining this error with the others. 

We next turn to the a('One blob' hadronic, 0{a^)) corrections, which are ob- 
tained by attaching a photon or fermion loop to the various lines in Fig. 1. They can be 



- 14 - 



THE HADRONIC CONTRIBUTIONS TO THE ANOMALOUS MAGNETIC MOMENT OF THE MUON 



further split into two pieces: the piece where both ends of the photon hne are attached to 
the hadron blob, o(h.v.p., 7), hadronic vacuum polarization corrections, shown in Fig. 6, 

and the rest. So we write, 

o('One blob' hadronic, 0{a^)) = a(h.v.p., 7) + a('One blob' hadronic, rest). (3.4) 

The last can be evaluated '^^1 in terms of the hadronic contributions to the photon vacuum 
polarization, finding 

10^^ X a('One blob' hadronic, rest) = -101 ± 6. (3.5) 

This result has been checked independently recently, by the Marseilles group (S. Priot, 
D. Greynat; E. de Rafael, private communication) and in ref. 21. 

The only contribution that requires further discussion is that depicted in Fig. 6, 
a(h.v.p., 7). In principle, this contribution can be evaluated straightforwardly by a gen- 
eralization of the method used to evaluate the O(a^) contributions. We can write 

/•oo 

o(2)(Hadr.) +a(h.v.p., 7) = / dt K{t)R^^\t), (3.6) 

i?(2)(t) 



where 

a^^\e'^e~ — > hadrons) +o''^'(e+e^ hadrons) + (J^°^(e"'"e~ — > hadrons; 7) 



CT(0)(e+e- -> At+A*-) 

The notation means that we evaluate the hadron annihilation cross section to second order 
in a, and we add to it the first order annihilation into hadrons plus a photon. For energy 
large enough this can be calculated with the parton model, which leads to a (very small) 
correction, (2 ± 1) x 10^^^. 

Then comes the contribution of small momenta. We start by discussing the process 
involving two pions. We calculate the corresponding piece by adding the contribution of 
the diagrams in Fig. 3, as given in the 2002 paper by Akhmetshin et al.I'*] In this way, we 
find 

10" X a(h.v.p., 7r+7r-7) = 47.6 ± 0.3. (3.7) 

The number is very close to that obtained in ref. 6 (46 ± 9) but the errors have decreased 
drastically.^ A similar analysis ought to be made, in principle, for other radiative interme- 
diate states like 37r-|-7 and KK + ^, which can be estimated in terms of the corresponding 
decays of the u) and </>, but they give a contribution below the 10^^^ level and we neglect 
them. The contribution from 7r''7r''7, (2.0 ± 0.3) x 10^^^, is taken from ref. 7. 

The lowest energy contributions to a^^\e'^e~ — > hadrons; 7) are those of the 
intermediate states 7r°7 and 777, Fig. 7. These contributions were evaluated in detail in 
ref. 6; we have 

10^1 xa(h.v.p., A) = 37±3 
10^^ xa(h.v.p., r?7) = 6±2 



® In ref. 6 only the radiation of hard photons weis included, evaluated using the results of ref. 19. 
The fact that this is so similar to the full result justifies the (expected) smallness of soft photon 
plus vertex correction. 
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Figure 7. The 7r°7, 

?77 contributions to 
o(h.v.p., 7). 



Collecting all of this, we get the total effect of the states hadrons + 7, including 
the loop correction, 

10^^ X a(h.v.p., 7) = 95 ± 6. (3.9) 

4. Byproducts 

4.1 The electromagnetic coupling on the Z, Q;q.e.d.(-^2) 

With a simple change of integration kernel the previous analysis can be extended to eval- 
uate the hadronic contribution to the QED running coupling, aQ.E.D.(^)j particular on 
the Z particle, t — M^; this is an important quantity that enters into precision evalua- 
tions of electroweak observables. By using a dispersion relation one can write this hadronic 
contribution at energy squared t as 



R{s) 



4„2 S(S - t) 

wlioro R is as in (2.1b) and the integral has to be understood as a principal part integral. 
Therefore, we can carry over all the work from the previous sections as in ref. 20, with the 
simple replacement 

K(s) — > — -. 7. 

We find, to next to leading order in a, 

10^ X Z\hada(M|) = 2 742 ± 12 (4.1) 
or, excluding the top quark contribution, 

10^ X Zihada^^^(Ml) = 2 749 ± 12. 
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Adding the known pure QED corrections, the running QED coupling, in the momentum 
scheme, is 

an EB (Ml) = . (4.2) 

y.ji.u.v z; 128.962 ±0.016 ^ ' 

The dilference with the result in ref. 20 is minute. (4.2) may be compared with a recent 
calculation of Hagiwara et al.,[^^l who, using only experimental e+e" annihilation data at 
low energy find 

aQ.E.D. (Ml) =^28.954 ±0,031- ^^"^^ 

While the central value is compatible with ours, the error in (4.3) is twice as large. This 
shows clearly the advantage of a thorough use of analyticity as well as of combining several 
sets of data (we will discuss this further in Sect. 6, in connection with the muon anomaly). 

4.2. The masses and widths of the rho, the quadratic radius and second 
coefficient of the pion, and the scattering length and effective range pa- 
rameter for the P wave in tvtt scattering 



As other byproducts of our analysis we can give very precise numbers for the masses and 
widths of the neutral and charged rho resonances, for the quadratic radius and second 
coefficient of the pion form factor, and for the the scattering length, a\, and effective 
range parameter, bl (defined as in ref. 9), for the P wave in nn scattering. In what regards 
the first, we have the following results: for the charged rho, and in MeV, 



TY, r, evr, e+e GJ (r decay only) 
mp± 774.0 ± 0.4 775.4 ± 0.4 

rp± 147.7 ±0.7 149.3 ±0.4 

X^/d.o.i. 1.06 ~ 1.8. 



(4.3) 



For the neutral rho. 



TY, e+e-,e7r TY, r, ctt, e+e" GJ (e+e" only) 

rripo 773.1 ±0.6 773.2 ± 0.4 772.95 ± 0.7 

Fpo 141.7 ±1.2 146.0 ±0.8 147.9 ± 0.7 

xVd.o.f. 0.9 1.06 ~ 1.7 



(4.4) 



Here TY are our results here, and GJ refers to the results by Ghozzi and Jegerlehner.P-^1 
Clearly, the masses are stable and well determined, the widths less so. 

The quadratic radius and second coefficient of the pion are defined by 

F2(^)^^^l + i(r^>^ + c,^^ (4.5) 

and our results allow a precise determination of both quantities. The same is true for a}, 
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6j. We have, for tt'^tt , 

TY [Only e+e^, ew] TY [Including r decay] Colangelo 

(rl) (W) 0.423 ±0.003 0.432 ± 0.001 0.435 ± 0.005 

c„iGeV^^) 3.78 ±0.05 3.84 ±0.02 

a\ {m^ = l) (38.9 ± 2.3) X 10-3 (37.8 ± 0.8) x lO^^ (37.9 ± 0.5) x lO^^ 
6} (m^ = 1) (4.1 ± 0.7) X 10-3 (4.74 ± 0.09) x lO^^ (5.67 ± 0.13) x lO^^. 

(4.6) 

TY is from our resiilts here and "Colangelo" from refs. 8, 9. The discrepancy between the 
result for b\ from the pion form factor and that in ref. 9 had already been noted in ref. 10. 
The numbers coming from our calculations for a\, bl, including tau decay data, agree much 
better than those not including it with results from pion-pion scattering, either from phase 
shifts analyses or using the Froissart-Gribov representation, I^*^! which is another reason 
for preferring the results including tau decays. 

The values of the parameters of our fits are 

e^e~, en : 

ci= 0.26 ±0.04, C2 = 0.19 ±0.13, (4.7) 
6o =1.106 ±0.009, 6i= 0.23 ±0.10 

and 

e'^e~ , CTT, r : 

ci = 0.24 ± 0.01, C2 = -0.18 ± 0.03, (4.8) 
6o =1.074 ±0.006, 6i = 0.13 ±0.04. 

The numbers hi in (4.8) correspond to tt^tt^; they are the ones we have used to calculate 
< 0.8 GeV^). For 7r°7r+ we would have, 

b(t = 1.064 ± 0.006, bl = 0.13 ± 0.03 [7r°7r+], (4.9) 

and the Ci are as in (4.8). The corresponding scattering length and effective range param- 
eter are, in units of m^+, 

a\ = (37.8 ± 0.8) x 10"^, bl = (4.78 ± 0.09) x 10"^ [tv\+]. (4.10) 



5. Comparison of our theoretical calculations with experiment 
for a„ 



We return to the magnetic moment of the muon, and present, in Table I, a summary of our 
results for a(Hadr.). In this Table we have added our old result from 2002, ref. 6, and the 
result of a recent evaluation by Jegerlehner,'^^] in which theory is kept to a minimum; for 
example, only data are used for s < 0.8 GeV'^, and QCD is taken valid only for energies 
above 13 GeV. The main interest of this type of calculation lies in its role as control of 
the calculations where a more comprehensive use of theory (as well as extra experimental 
information) is made. In the Table we also include the recent experimental value of ref. 2. 
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For our evaluations here, we have separated expUcitly the estimated errors due to radiative 
corrections, and the hght-by-hght scattering diagram. 







10" X ((^(Hadr.) 


Jegerlehner 




6 840 ±94 


TY (2002), e+e", eTr, t 


4 774 ±51 


6 993 ±69 


TY, e+e", ctt 


4 715±32± 10(rad.) 


6 935 ± 50 ± 10 (rad.) ± 30 (£ x £) 


TY, e+e", eTr, r 


4 798 ±31 ±10 (rad.) 


7018 ± 49 ± 10 (rad.) ± 30 (^ x I) 


Experiment 




7 209 ±60 



Table I 

Contributions to the rho region, and to the hadronic part of the muon anomaly. Jegerlehner: 
ref. 22. TY (2002): ref. 6. TY: this article. 



Erom Table I it is clear that there is reasonable agreement among the various 
theoretical determinations, but there is a definite distance between the central values 
from theory and experiment, at a level between 2.3 ct and 3.3 a, if we add quadratically 
"rad" and "£ x £" errors to the other ones. In the remaining of this section we will discuss 
possible reasons for this discrepancy. 

An obvious reason would be new physics; we will not discuss this here, since it 
lies outside the scope of the present paper, and send the interested reader to the hundreds 
of papers that have been written discussing this possibility. 

A second reason is, of course, a displacement of the experimental result. Since 
the experimental number in Eq. (1.1) comes basically from only one experiment, it could 
happen that an independent determination would move it to better agreement with the 
results of the theoretical evaluations. 

And a third possibility is that the central values of some of the theoretical eval- 
uations presented here are displaced with respect to the true values. The more obvious 
place where such a displacement may occur is the evaluation of a(^ x I). The two approx- 
imations used to evaluate this contribution do not have overlapping ranges of validity; 
there is a wide region, when the virtualities of some (one) of the intermediate photons 
are small, and at the same time other (others) are large, where neither the one-pion or 
the constituent quark approxiuiations need to be valid. In fact, the only certain result 
we have on this piece is the coefficient of the leading chiral logarithm (log'^ m^) given 
in ref. 23, which is of little practical use. However, it is not easy to see how one could 
get the large values necessary for theory and experiment to overlap. In a recent calcu- 
lation, using methods somewhat different to previous ones, Melnikov and Vainshteinl'^^1 
find 10^^ X a(£ x £) = 136 ± 25. Although this is larger by a bit more than one sigma 
(44 X 10~^^) than the result quoted in (3.3), it is not sufficient to remove the discrepancy. 
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Finally, we may have a coincidence of several of the effects mentioned here, with 
the bad luck that they add. 

6. Comparison with other recent calculations and 
concluding remarks 

Our aiialj'sis shows that, to get a precise value of a^, it is certainly necessary to profit 
from the existence of methods that allow us to make full use of theory, in i)articular for 
fitting the pion form factor: the use of robust theory in the fits prodiices robust results. 
It is not very consistent to use analyticity and unitarity to write the representation (2.1), 
and refuse to use exactly the same ingredients to improve the knowledge of F^r. In this 
sense, it is also important to take into account the data on Fj^{s) for spacelike s, i.e., 
from en scattering. This has been deemed inappropriate by some authors because they 
(may) contain systematic errors. However, if one refrained from taking into account data 
affiicted by systematic errors, there would be no data one could use. Systematic errors 
can and should be taken into account as shown, in this particular instance, in ref. 6 and 
in the present article. 

In the same vein, we believe that tau decay data can and should be used, in spite of 
the fact that the data of the various tau decay experiments differ in some energy regions 
by more than one standard deviation -doiibtlessly because of systematic errors, as is 
obvious from (2.11c). This is particularly important because the errors given in Eq. (2.8a) 
for a'^\s < 0.8 GeV^) are deceptively small. As we already commented, the error per 
experimental point for e+e~ — > tttt, fitting only e+e^, en data, is 89/113 ~ 0.79, clearly 
smaller than unity; while, even imposing tau decay information, the error per experimental 
point is (as reported in Eq. (2.11c)) only of 108/113 ^ 0.97, perfectly acceptable. This 
means that acceptable fits -like, indeed, the one obtained by us using also tau decay data- 
can be found outside the nominal error bars in Eq. (2.8a).® 

Likewise, when data on from the processes e+e^ tt+tt^ +7 are forthcoming, 
they should be incorporated into the analysis: the safest way to get rid of systematic errors 
is to combine data of various, independent experiments, so that the various independent 
systematic errors average out. The gain, both in accuracy and robustness that follows 
from our methods can perhaps be seen more clearly if we compare them with other recent 
evaluations, P^'-^^'^^l something that we do in Table II (where we do not include the results 
of ref. 22, already discussed before). Here the stability of our 2002 results against including 
new, more precise e+e" data, contrasts with the variations in the other determinations, 
is spite of the fact that our error is substantially smaller. 

In short: fitting to the theoretical expressions instead of integrating directly the 
data, allows us to compare the different data samples among them in a fully quantitative 



One should, however, not forget that the value we needed for the tau decay normalization error, 

(-1.4±0.5)% (cf. Eq. (2.11b)), or (-1.2±0.5)% if including the estimate (2.16) for 5^, is slightly 
larger than the expected normalization error, as given by the Particle Data Tables, 0.7%. 



-20- 



THE HADRONIC CONTRIBUTIONS TO THE ANOMALOUS MAGNETIC MOMENT OF THE MUON 





10" X a^""^' ^°+^ 


10" X 


Ezhela [e"''e~] 


6 996 ± 89 


116 591 835 ±96 


Hagiwara [e'''e~] 


6 924 ± 64 


116 591 763 ±74 


Davier (a) [e"'"e~] 


6 847 ± 70 


116 591 693 ±78 


Davier (b) [e~'"e~] 


6 963 ± 72 


116 591 809 ±80 


Davier (b) [r] 


7 110 ±58 


116 591 956 ±68 


TY (2002), e+e", ctt, r 


7 002 ±66 


116 591 849 ±69 


TY, e+e", ctt, t 


7 027 ±49 


116 591 889 ±58 


Experiment 




116 592 080 ±60 



Table II 

The lowest order hadronic part of the muon anomaly, including photon corrections [Eq. (3.9)], 
and the full ai_i. Hagiwara: ref. 21. Ezhela: ref 25. Davier, (a) and (b): ref. 26, (a) and (b). TY 
(2002): ref. 6. TY: this article. 



manner. We believe that only using this kind of quantitative comparisons one can decide if 
the suggested discrepancies are meaningful or only apparent, and, if meaningful, if they are 
due to systematic errors or to physics. The fact is that one can fit all the en, e+e^ mr 
and tau data with a x^/d.o.f. essentially 1. This result is a non-trivial improvement on 
previous work. 

It is, however, not clear to us that one can improve the results using mr scattering 
data. If we include in the fit the experimental numbers for 5{ , the value of increases by 
8 X 10~^^; but this is not necessarily more precise than the result without including this 
information. The systematic errors of 5l, due to the fact that one does not scatter real 
pions (and thus one has to rely on models), are larger than the errors in our calculation. 
And if, like Colangelo and collaborators, t^'^l we input 6} from theoretical analyses (Roy 
equations and chiral perturbation theory), one is depending on determinations whose 
accuracy has been challenged!^"! and is, very likely, too optimistically estimated. 

Apart from this, to get real improvement in the theoretical predictions for the 
quantity a^(Hadr.) it would be, first of all, necessary to remove the sources of uncertainty 
mentioned at the end of Sect. 5. Of these, the one stemming from £ x £ is unlikely to be 
removed in a satisfactory manner; and we have also a problem (although less important 
numerically) with electroweak radiative corrections to tau decay. If we treat them by 
considering the pions as elementary, and factoring out F^, then Sirlin's theoreml^^l implies 
that, for T umr, the logarithmic piece, log Mz/rrij- cancels out. It, however, does not 
cancel if we consider that, at short distances, the decay is really t~ — > I'iid. This means 
that the model with elementary pions fails, for this case of tau decay, and thus that the 
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chances of removing the uncertainties due to lack of accurate knowledge of the radiative 
corrections here are remote. 

As a final comment, we would like again to bring attention to the mismatch 
between the experimental and theoretical values for a^; although not yet definite proof of 
failure of the standard model, it cannot be ignored. 
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